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We analyze the phemomenon of induced fermion number at finite temperature. At finite temper- 
ature, the induced fermion number (A^) is a thermal expectation value, and we compute the finite 
temperature fluctuations, (AA'^)^ = (A^'^) — {N)^ . While the zero temperature induced fermion 
number is topological and is a sharp observable, the finite temperature induced fermion number is 
generically nontopological, and is not a sharp observable. The fiuctuations are due to the mixing 
of states inherent in any finite temperature expectation value. We analyze in detail two different 
cases in 1 + 1 dimensional field theory: fermions in a kink background, and fermions in a chiral 
sigma model background. At zero temperature the induced fermion numbers for these two cases 
are very similar, but at finite temperature they are very different. The sigma model case is generic 
and the induced fermion number is nontopological, but the kink case is special and the fermion 
number is topological, even at finite temperature. There is a simple physical interpretation of all 
these results in terms of the spectrum of the fermions in the relevant background, and many of the 
results generalize to higher dimensional models. 



I. INTRODUCTION 



The phenomenon of induced fermion number arises due to the interaction of fermions with nontrivial topological 
backgrounds (e.g., solitons, vortices, monopoles, skyrmions), and has many applications ranging from polymer physics 
to particle physics ||| |, |, |. The original fractional fermion number result of Jackiw and Rebbi [Q provides the 
physical explanation of the existence of spinless charged excitations in polymers ||^, ^, |^, ^, ^ . The adiabatic analysis 
of Goldstone and Wilczek in systems without conjugation symmetry has important implications for model field 
theories in particle physics, such as bag models [^, electroweak theories and chiral sigma models [^. At zero 
temperature, the induced fermion number is a topological quantity, and is related to the spectral asymmetry of the 
relevant Dirac operator, which counts the difference between the number of positive and negative energy states in 
the fermion spectrum || | , |, |l3|, [ll , [Tsl , [T7| , p^ , [l9[ |. Mathematical results, such as index theorems |Q, ^ and 
Levinson's theorem ||l5, 17|, |18[, imply that the zero temperature induced fermion number is determined by the 
asymptotic topological properties of the background fields. This top olog ical character of the induced fermion number 
is a key feature of its application in certain model field theories |10, pT[ p2[. 

At finite t emp erature, the situation is less clear-cut. In several examples, namely the 1 + 1 dimensional chiral kink 
background |po[[2l], and the 3 + 1 dimensional Dirac ^ and 't Hooft-Polyakov monopolc |^ backgrounds, the 
finite temperature induced fermion number has been shown to be temperature dependent, but still topological in the 
sense that the only dependence on the background field is through its asymptotic properties. However, it has recently 
been demonstrated [ p6[ that the finite temperature induced fermion number need not be topological. Explicitly, in a 
' 1 + 1 dimensional chiral sigma model, the finite temperature induced fermion number depends on the detailed structure 
of the background, not just its asymptotic behavior. This result corrects several previous analyses [ p7| p8| that had 
claimed that the finite temperature fermion number was in general a topological quantity (at zero chemical potential). 
In [^6| , an explicit calculation showed that the finite temperature fermion number can be nontopological, and a simple 
physical explanation was given for the origin of the nontopological fermion number as the plasma response of the 
fermions to the inhomogeneous background. A possible source of confusion here is that while the kink and sigma 
model cases are very similar at zero temperature, it is not widely appreciated that at finite temperature they are very 
different. In this paper we present more details of this analysis, and we present a general argument that the finite 
temperature fermion number naturally separates into a temperature-independent topological piece that corresponds 
to vacuum polarization effects, and a temperature-dependent piece that is generically nontopological, and which 
corresponds to the thermal occupation of excited states of the fermion Fock space. While this temperature-dependent 
piece is generically nontopological, for certain particular backgrounds (such as, for example, the kink background) 
the fermion spectrum has a special symmetry, which has the consequence that the temperature-dependent corrections 
are in fact themselves topological. 

Another motivation for our work concerns the nature of the induced fermion number as a sharp quantum observable. 
At zero ternperature it has been shown that the fractional induced fermion number is indeed a sharp observable 
p9| , |30| , [si] , p2| , |33| . In this paper we address this question at finite temperature by computing the rms fluctuations, 
{ANf — {N^ — (N)^, in the induced fermion number (N). We compute the finite temperature fluctuation in two 
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different cases, and show that AA'' is generally nonzero (and nontopological) at finite temperature, but that AN 
vanishes at zero temperature. This indicates that N is not a sharp quantum observable at nonzero temperature. 
The nonvanishing fluctuations are due to the mixing of Fock states inherent in the thermal expectation value. It is 
interesting to note that another example of non-sharp fractional fermion number (but not in the context of temperature 
dependence) has been discussed recently in the context of liquid helium bubbles |Q . 

Our analysis of finite temperature fermion number has also been motivated by the realization over recent years 
that certain types of anomalies in zero temperature field theory become much more subtle at finite temperature. For 
example, Pisarski et al have explained ||35|| why finite temperature anomalous 7r° decay amplitudes are temperature 
dependent even though the chiral anomaly (a topological object that is known to be related to the anomalous 7r° 
decay amplitude at T = 0) is known to be temperature independent. These and related issues have also been explored 
in 1 + 1 dimensions |3^. And in odd spacetime dimensions, it has recently been realized |3^, ^ Q that while 
the topological Chern-Simons term is the only parity-violating term that can be induced in the T — effective action, 
at finite temperature there are infinitely many other parity-violating terms, all of which vanish identically at zero 
temperature, but all of which are crucial, for example, for understanding how it is possible to maintain large gauge 
invariance at finite temperature. These results have led us to reconsider the related general question of induced 
fermion number at finite temperature. 

In Section II we define what is meant by finite temperature induced fermion number, and indicate how it can be 
computed. In Section III we give more details of the derivations and results of p6[ for the finite temperature induced 
fermion number (N) for a kink background and for a sigma model background in 1 + 1 dimensional field theory. 
Section IV presents the computation of the finite temperature fluctuations, (AiV)^ — {N'^) — {N)'^, in the induced 
fermion number for the kink case and for the sigma model case. In Section V we give a simple physical interpretation 
of all these results in terms of the thermal occupation of the available levels in the Dirac spectrum, according to 
Fermi-Dirac statistics. We conclude in Section VI and give some general comments regarding the extension of our 
results to higher dimensional theories. 



II. FINITE TEMPERATURE INDUCED FERMION NUMBER 



Consider an abelian model in 1 + 1 dimensions with fermions interacting via scalar and pseudoscalar couplings 
to two bosonic flelds i/ii and 02. For the purposes of this paper, and (j)2 will be considered as static classical 
background fields. The Lagrangian is 

C = i-tpp'ip-'4){<pi+ij5(l)2)'ip (1) 

We will concentrate on two important physical cases: 

(i) kink case [Q: 0i is constant, but 02 = <l>2{x) has a kink-like shape, as shown in Fig. 1: 

01 = m, and 02 = kink, with 02(±cxd) — ±02 (2) 

(ii) sigma model case ||^: 0i and 02 are both x-dependent, but are constrained to the "chiral circle": 

<Pl+<Pl= (3) 
In each case, we define a corresponding angular field 

e{x) = arctan ( (4) 

In the kink case (0), the angular field 9{x) = arctan (^|^) also has a kink-like shape, as shown in Fig. 1. In the 
sigma model case (0), the angular field 9{x) has the interpretation of a local chiral angle, since the chiral constraint 
(P|) allows us to write (pi — m cos 9 and (j)2 — ra sin 0, so that the interaction term in the Lagrangian (Q) is 

m^/i (cos 9 + sin 9)ip = e^^'^^ tp (5) 

In this sigma model case, it is the angular field 9{x) that has a kink-like shape, as shown in Fig. 2. 

The second-quantized fermion number operator is defined |^ as = ^ , VP] . The fermion field operator ^ 
may be expanded in a complete set of eigenstates of the Dirac Hamiltonian H for fermions in the presence of the 
background fields 0i and 02. The Dirac Hamiltonian corresponding to the Lagrangian (|l|) is 

H = -i7°7^V +7"0i(a;) +«7*'7502(a;) 

01 -iV-i02\ 
-^V + 02 -01 / ^' 
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FIG. 1: In the kink case l)^, the kink function 4>2{x) has the form shown in the sohd line, with asymptotic values ±(j)2- The 
corresponding angular field defined in (^), 9{x) = arctan(^^^), also has a kink-like shape, as shown in the dashed line, with 

asymptotic values pmO, where 9 — arctan(^). Note that in the kink case, the allowed values of 9 lie in the first branch [—-J, ^] 
of the arctan function. 

where V = and we have chosen to work with the Dirac matrices: 7" = (T3, 7^ — ia2, and 7^ — — (Ti. The 
presence of the background fields modifies the fermion spectrum from the free case. All necessary information about 
the fermion spectrum can be encoded in terms of the spectral function (j{E) of the Dirac Hamiltonian in (0): 



a(E) = - Im Tr f 

The zero temperature vacuum expectation value, (0|A^|0) = (A^)o, of the number operator can be expressed as 



dEa{E) sign(i;) 



(7) 



(8) 



where the subscript refers to zero temperature. Thus, {N)o is expressed in terms of the spectral asymmetry of the 
Dirac operator, which counts the number of positive energy states minus the number of negative energy states. 




FIG. 2: In the sigma model case (^), it is the angular field 9{x) defined in (^) that has a kink- like shape, with asymptotic 
values ±9. 
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At zero temperature, in both the kink case ^ and the sigma model case (^), the induced fermion number (||) is 
given by the simple expression |^ : 

(iV)o = — / dx0'^-e (9) 

where 9{x) is the angular field defined in (Q), and 9 is the asymptotic value of 9{x) at x = +oo. The zero temperature 
fermion number (A^)o is topological in the sense that it depends only on the asyinptotic value 0, and not on the detailed 
shape of 9(x). The original conjugation symmetric case of Jackiw and Rebbi W is obtained by taking m ^ in the 
kink case (||), in which case (A^)o ^ i^- Section IV we will show that the zero temperature expectation value (H) 
is in fact also a sharp eigenvalue. 

At nonzero temperature, the induced fermion number is a thermal expectation value 



{N)t - 



TT{Ne-P") 
Tr {e-(>") 

2.1 ^ ' ' V 2 



^ ^ dEcr{E) tanh ] (10) 



where P — is the inverse temperature. Note that as T ^ 0, which means — *■ oo, only the vacuum state survives 
in the trace, and {N)t reduces to the vacuum expectation value (iV)o- Correspondingly, tanh(^) sign(i<^), so 
that the integral expression in ( p^ reduces smoothly to the spectral asymmetry expression in (||). In fact, the finite 
temperature expression ( [l0| ) provides a physically natural, and computationally simple, regularization of the spectral 
asymmetry (^. In Section IV we will show that at nonzero temperature the fermion number expectation value ( p^ ) 
is not a sharp eigenvalue. 

To compute {N)o or {N)t, one needs information about the spectral function cr{E). We stress, of course, that (t{E) 
has nothing to do with the temperature; it simply describes the spectrum of the fcrmions in the presence of the static 
background fields 4>i and <^2- One convenient way to proceed is to use the expression (Q) for the spectral function to 
write ( px| ) as a contour integral in the complex energy plane: 

Here tr(-^ji^) is the resolvent of the Dirac Hamiltonian H, and the contour C is (— oo + ze, +oo + ie) and (+oo — 
ie, —oo — ie), as shown in Fig. 3. Since H has a real spectrum, we can evaluate the contour integral (|ll| ) in two 
alternative ways. First, we can deform the contour around the simple poles of the tanh(^) function, which occur 
along the imaginary axis at the Matsubara modes: z„ = (2n + l)i7rT, for n £ 7i. This leads to an expression for 
{N)t as an infinite sum. Alternatively, we can deform the z contour around the poles and cuts of the spectrum of H, 
which lie on the real axis. This leads to an integral representation for {N)t, which is just the Sommerfeld- Watson 
transform of the infinite sum expression. We will see examples of these equivalent forms below. 



III. RESULTS FOR FINITE TEMPERATURE INDUCED FERMION NUMBER 



Since tanh('^) is an odd function of E, it follows from ( ]lO| ) that to compute {N)t we only need the odd part, 
o'odd{E) = ^[(j(E) — (t{—E)], of the spectral function: 



/•oo 

(iV>T - - / dE<JoM{E) tanh 

"'0 



f) 



(12) 



Correspondingly, we only need to know the even part of the resolvent: 



(^)t = - 



dz 
27ri 



tr 



1 



H - z 



tanh 



Pz 

T 



(13) 



This fact simplifies the calculations considerably, and has important physical consequences, as will be discussed in 
Section V. 
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FIG. 3; The contour C in the complex energy plane used in the evaluation of the contour integral in ([ill). The spectrum of the 
Dirac Hamiltonian lies on the real axis and is indicated schematically by the continuum cuts and some possible bound states 
in the gap. The crosses on the imaginary axis depict the simple poles, z„ — (2n + l)7rir, of tanh(^). 
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FIG. 4: Plots of (tt times) the finite temperature induced fermion number {N\ given in ( |l7| ) and (^_8|), for the kink case. The 
plots show the dependence of 7r(A'^) on the asymptotic value B = arctan(^), for various values of the temperature, as indicated. 
Note that for low temperature (large (A'^) smoothly approaches the zero temperature result (A'^) = ^/tt in (^. 



Finite T Fermion Number in the Kink Case 



For the kink case (||) we can use the remarkable result (a special case of the Callias index theorem [^) that the 
even part of the resolvent (and hence the odd part of the spectral function) can be computed exactly for any kink 
02 (x): 



tr 



R-z 



-ni(f)2 



even (^2 _ ^2) 



(14) 



Notice that this only depends on the kink 4'2{x) through its asymptotic value 02! it is otherwise independent of the 
details of the shape of the kink. 
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FIG. 5: Plots of (tt times) the finite temperature induced fermion number (A'^), given in ( |l7| ) and (y_8|), for the kink case. The 
plots show the dependence of 7r(A^) on the inverse temperature for various values of the asymptotic value Q, as indicated. 

Note that for low temperature (large /3), {N) quickly saturates to the zero temperature value 7r(A^) = d. Also note that (A'') 
vanishes as T cx). 



To motivate this fact, we can express the Dirac resolvent as tr (^75^ j = tr yji + z)jj^—^^ , and note that, for the 
kink case, iJ^ takes the simple diagonal form 



if2 = 



Thus, 



tr 



H 



tr 



- tr 



(15) 



(16) 



That is, we can express the even part of the Dirac resolvent as the difference of two Schrodinger resolvents, with 
the corresponding Schrodinger potentials, V± = (/)2 ± 02, being isospectral This isospectrality property is the key 
ingredient for proving (|l^ . 

Given the result (^^, the induced fermion number (|l^) for the kink case (||) is 



=0 ((2n + 1)2 + (1^)2) y^(2n + l)2cos2^+(^)2 



(17) 



where we recall that 9 = arctan(^), and we can restrict our attention to < < ^. As an alternative to the 
summation expression in ([l7|), we can write (A'^)t as an integral representation: 



1 . , yiTid. sin 6* cos f 
(7V)t = - sign(0) tanh(-^) 

Z Z TT 



du ■ 



tanh(^^f sec6'u) 



\At2— T {v? — C0s2 6) 



(18) 



The kink case induced fermion number (iV)T is plotted in Fig. 4 as a function of 9 for various values of the temperature, 
and in Fig. 5 as a function of for various values of 9. The integral form (|lj) makes the zero temperature limit 
clear : using the integral 



du ■ 



\/u2 — 1 (v? — cos2 



sm W cos ( 



(19) 



it follows that the leading correction to (iV)o at low temperature is exponentially small : 



sign(6') e 



-7nf3 



f3^ 00 



(20) 
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Thus (-/V")t reduces smoothly to the zero temperature expression (^. This can be seen clearly in Figs. 4 and 5. But 
at finite temperature, the expressions ( |l7| ) and (18) for the induced fermion number are much more complicated than 
the T = result (^. Nevertheless, (-/V)t in ( p^ ) and ( p^ is still topological in the sense that it only depends on the 
kink background through the asymptotic value 6 — arctan((/)2/m). Other details of the kink shape do not matter. 
We will see in the next subsection that this is not true in the sigma model case (|^). It is also interesting to note that 
the angular nature of the parameter 9 is clearly manifest in the finite temperature expressions ( p7| ) and (p^), while 
it is not as obvious from looking at the zero temperature expression (^. A similar observation applies for the finite 
temperature induced fermion number in 3 + 1 dimensions for fermions in the presence of a Dirac monopolc |2^, ^ , 
and for fermions in the presence of a static 't Hooft-Polyakov monopole [ESl Em . 



B. Finite T Fermion Number in the Sigma Model Case 

In the sigma model case (||), the Callias index theorem result (|lj) for the even part of the resolvent does not apply. 
There is no general expression for the even part of the resolvent. This is because in the sigma model case 

i/2 = f-^' + "' 2]+me' i':^'^^ -'''''I] (21) 

\^ —V +m^j \ismd —cos6 J ^ ' 

which is clearly not of the diagonal isospectral form in (|l5|). Thus, another approach is needed to evaluate the 
resolvent. In |26[ , the derivative expansion was used to evaluate the even part of the resolvent as an expansion in 
powers of & and its derivatives. The derivative expansion Q assumes that 0' <^m : the spatial derivatives of the 
background fields are assumed small compared to the fermion mass scale m. In other words, the background chiral 
field 0{x) is assumed to be slowly varying on the scale of the fermion Compton wavelength. To next-to-leading order 
in the derivative expansion for the sigma model case, one finds [ p6[ 



tr 



1 



H 



2(m2 - z2)3/2 J 8(m2 - z2)5/2 J I6(m2 - z'^yi' 



where the dots refer to terms involving five or more derivatives. The middle term vanishes since J dx9'" = if 9(x) 
has a kink-like shape. Inserting this approximate expression into (|l^) we obtain: 

WT = -(-j IS [(2„ + 1,^ + (H.,.p/. / S |(2„ + 1). + y 

The first term is topological and reduces smoothly to the zero temperature result (^ as T — > 0. The second term 
involves / dx {9'Y , which is clearly not topological : it depends on the actual shape of the chiral field 0{x), not just its 
asymptotic value 9. This is still consistent (to this order) with the topological nature of the zero temperature induced 
fermion number (|^), because the energy trace prefactor multiplying the J dx {9'Y term in ( [2^ ) vanishes at T = 0. 

Higher orders in the derivative expansion ( |2^ ) can be developed systematically, although it becomes somewhat 
tedious to enumerate all the different terms at high orders. However, a remarkable simplification occurs [p6| in the 
low temperature limit where T <C m. The leading low T terms at each order of the derivative expansion have the 
simple form : 



21-1 



Thus, in the low temperature limit, we can resum the entire derivative expansion, to obtain the induced fermion 
number in the sigma model case (0) : 



where the dots refer to power-law subleading terms for T <^ m. 

Note that {N)t in ( |25|) naturally splits into a temperature independent topological term, and a temperature 
dependent nontopological term. This has a simple physical explanatio n |B6p . First, observe that the chiral sigma 
model background acts like a spatially inhomogeneous electric field |l^, [l4[, as can be seen by making a local chiral 
rotation: ^ ^ ijj = In terms of these chirally rotated fields the Lagrangian (|l|), with interaction (||), becomes 

- _ - _ - 9' ~ 
£ = iip^ — mip'4' ~ 7° "2" V' (26) 



8 



Thus, the chiral field 9{x) acts as an inhomogeneous scalar potential Aq{x) = ^9'{x), leading to an inhomogeneous 
electric field 

Eix)^-^e'\x) (27) 

This electric field acts on the Dirac sea to polarize the vacuum by aligning the virtual vacuum dipoles of the Dirac 
sea, producing a localized build-up of charge near the kink center. This leads to the first, topological, term in ([2^), 
which is just the familiar zero temperature result [^3[ It is temperature independent as the short-lived virtual 
electron-positron dipoles of the Dirac sea do not come to thermal equilibrium. The second, nontopological, term in 
( p5| ) arises as the response of the real charges in the thermal plasma to the spatially inhomogeneous electric field . 
The linear response ||44| of the plasma at low temperature to such an electric field yields an induced fermion number 
density 



/lik 
(28) 



where f{x,k) = f+{x,k) — J-{x,k) is expressed in terms of the particle and antiparticle Fermi densities. In the 
derivative expansion limit, the leading effect of the background is that of an inhomogeneous chemical potential 
IJ.{x) = — ^ 0' [x). Thus, in this limit, the local Fermi particle and antiparticle distribution functions are 

" e/3(yfe5WTP(x)rrT ^^^^ 

The resummed derivative expansion expression (p5| ) was obtained in the low temperature (T <^ m) limit. In this 
limit, we can write 



f+{x,k) - Ux,k) « ^-PV^W ^ ^-PK^)-^ 



= -2e-^vW j^^ j (30) 

Note that f+ — /_ vanishes smoothly as T ^ in the derivative expansion regime, because 9' <ti m. For low 
temperature, the leading behavior of the k integral in (psh is 



^e-P^^^ ^ J—e-^"^ + ... , T^O (31) 
, 27r V 27r ^ ^ 

Therefore, the plasma linear response contribution ( p8| ) to the induced fermion number gives precisely the second, 
nontopological, term in the formula ( p5| ) which was obtained by resumming the derivative expansion at low temper- 
ature. 



IV. THERMAL FLUCTUATIONS OF INDUCED FERMION NUMBER 

We now turn to the fluctuation, (AA^)|, = {N'^)t — (^)ti ^^'^ induced fermion number {N)t- We will show that 
the fluctuation vanishes at zero temperature, but is nonzero at flnite temperature. Furthermore, the fluctuation is 
inherently nontopological, in both the kink and sigma model cases. 

Recall that the fluctuation can be expressed in terms of the partition function as |^ 

Thus, the fluctuation can be expressed in a manner analogous to the fermion number expressions (^0|) and ([III): 

2 ' 



3sed in a manner analogous to the 
(AN)^ = - / dE(j{E)sech^{- 

4 J —oo 



I f dz ^ f 1 \ ^ 2 f 



■ tr sech^ — (33) 



From a purely computational viewpoint, there is an immediate difference between ( p3[ ) and the fermion number 
expressions ( p^ ) and ( pT| ) : since sech^(^) is an even function, we now need the even part of the spectral function 
(equivalently, the odd part of the resolvent), whereas to compute (-/V)t one needs the odd part of the spectral function 
(equivalently, the even part of the resolvent) . 
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A. Fluctuations in the Kink Case 



The Callias index theorem resuh (|Tj) gives the exact form of the even part of the resolvent, but tells us nothing 
about the odd part of the resolvent. Thus, in the kink case, we should no longer expect an exact topological result for 
the fluctuation (AN)'^. Indeed, for a general kink background, an approximation is needed to compute the odd part 
of the resolvent. Instead, here we will proceed by considering a special two-parameter family of kink backgrounds for 
which the both the even and odd parts of the resolvent can be computed exactly. Consider 



02 (a^) = tanh 




(34) 



where j = 1, 2, 3, . . . is an integer. The parameter represents the asymptotic value of the kink, while j determines 
the scale of the kink, as illustrated in Fig 6. While not completely general, the family of kinks in ( |3^ ) is sufficiently 
general to distinguish between the topological and nontopological effects. For these special kinks, the isospectral 
potentials appearing in the square ( |l5| ) of the Dirac hamiltonian are 



± 



l^sech^^; 



(35) 



which are exactly solvable reflectionless Poschl- Teller potentials, for which the resolvents are known in closed form, 
as is most easily derived from the exact phase shifts [^6[ |48j. In fact, for the family of kink fields 02 (a;) in (|34|), 
the associated Schrodinger resolvents are : 



tr 



tr 



E 



^Jm^ + 4>^-z^ 1=1 {rn'^ + 0^(1 " jl) " z^) 



E 



i/j 



z- 1=1 (m^ + r(l-^)-z' 



(36) 
(37) 



Note that the only difference between (36) and ( pT]) is the I — j term in the sum, which is present in ( |37| ) but not 
in (pq). Since the even part of the Dirac resolvent (16) involves the difference of these two isospectral Schrodinger 



resolvents, we see that this difference is indeed independent of the scale j, and confirms the topological Callias index 
theorem result ( |l^ for this case. On the other hand, the odd part of the Dirac resolvent depends on both the 




FIG. 6: Plots of the special kink backgrounds in (jsj) for various values of j, illustrating the interpretation of the parameter j 
as a scale parameter that characterizes a nontopological aspect of these kinks. 
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FIG. 7: Plots of the fluctuation AA'^, given in (M) and ([40|), for the special kink backgrounds in (|34j). The plots show the 
dependence of AA^ on the inverse temperature ^^^~for various values of the kink scale j, as indicated. This plot is for 6 — tv/S. 
The solid line shows the corresponding value of the fermion number (N) , which is topological and therefore independent of the 
scale j. Notice that AiV vanishes exponentially as the temperature approaches zero, while (N) saturates to its zero temperature 
value of 1/8. These plots show the dependence of the fluctuation AA'^ on the scale j, indicating that AN is nontopological. 



asymptotic value (f) and the scale j: 
1 



tr 



H~ z 



\<t>\z 



Ddd 



i-1 



l/j 



(38) 



The modulus appears in the Dirac resolvent ( 38h because the complex Dirac energy z is the square root of the 
Schrodinger energy parameter z'^ appearing in (^) and (^), and so we must be careful about the sheet structure in 
the complex z plane. For example, from the Dirac Hamiltonian (p) it follows that if (/) > then there is an unpaired 



bound state ip — 







exp(- /(/)2) 



aX E = — m, while if < then there is an unpaired bound state "0 



_ / exp(/02) 







at E = +m. This fact is encapsulated in the Dirac resolvent (|lj) and (|38|), and similar arguments apply for the other 
bound states and for the continuum cuts. 

Given this result (^8|) for the odd part of the resolvent, the contour integral expression ( p3| ) for the fluctuation 
yields a summation expression: 



(AiV)^ 



2ml3\sin§\ 



E 

n=0 



2{2n+l)^cosH-{^)^ (!^)2_(2n + l) 



[(^^)2 + (2n+l)2] [(^i^)2 + (2n+l)2cos2 



3/2 



^i^n\ 2(2n + 1)4 cos2 0- (1^)2 (1^)2 _(2n+ 1)2 _ (IIlH)4(i _ tan^ 



(^)2 { 1 + (1 - 4) tan2 9) + (2n + 1)2 (1^)2 + (2n + 1)2 cos2 6 



3/2 



(39) 



Alternatively, there is an equivalent integral representation expression: 



(AiV)^ = ^ <; sech^(^) 



21 sin ^ 



TT 



u sech' 



du 



\ 2 



sec9u 



+2E 

1=1 



sech 



m(3 



'l + (1 - p)tan2 



2|sin6l| 

J 



TT 



(u2 — C0S2 ^)\J — \ 

usech^ (^lil^sec^ii) 



du 



(^.2_( 



1 — 4r sm 



^ .in^ 9) sf^F^ 



(40) 



The fluctuations are plotted in Fig. 7 for = ^, and various values of the scale parameter j. It is clear that the fluc- 
tuation vanishes exponentially at zero temperature (m/3 co). It is also clear that the fluctuation is nontopological. 
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as it depends on the scale parameter j. At low temperature, the leading behavior is independent of the scale j 

(A7V)2, ~ e^""^ + ..., (41) 

but the subleading exponential corrections are j dependent, as can be seen from the integral expression ( jiol ) and from 
the plots in Fig. 7. 



B. Fluctuations in Sigma Model Case 

In the sigma model case, there is no exact expression for either the odd or the even part of the Dirac resolvent. 
However, applying the derivative expansion as in [ p6| , we find that in the low temperature {T <C m) limit the dominant 
contribution comes from terms involving (even) powers of 9' [x). These can be evaluated and resummed to all orders 
of the derivative expansion, to yield the leading low temperature result : 



[AN f^ ^ e-"/^ j dx sinh^ (^) + ■ • ■ (42) 

This expression has a smooth low temperature limit in the derivative expansion regime where 9' <^ ni. Indeed, the 



fluctuation (42) vanishes as T — > 0. Furthermore, at any finite temperature, the fluctuation ( [42[ ) is clearly nontopo- 
logical. In the next Section we will see that this expression for the fluctuation has a simple plasma interpretation, 
analogous to the linear response explanation (|2|) - (|3l| ) of the induced fermion number ( |25[ ) in the sigma model case. 



V. PHYSICAL INTERPRETATION 



In this Section we give a simple physical interpretation of the results of the previous Sections for the induced 
fermion number and its fluctuations, in terms of the spectrum of the fermions in the background field. Let us first 
summarize the results. At zero temperature, the induced fermion number (A^)o is always topological, and has vanishing 
fluctuation: (AiV)o = 0. At finite temperature, (7V)t is topological for a kink background, but nontopological for a 
sigma model background. For both the kink and sigma model cases, the fluctuation is non- vanishing, (AA'^)t 7^ 0, 
at nonzero temperature, and (AA'^)t is nontopological. The vanishing of the fluctuations at zero temperature is 
in agreement with previous arguments that the zero temperature induced fermion number is a sharp observable 
|p9| , po| , ^ 33 . At finite temperature, the nonvanishing fluctuation indicates that the induced fermion number is 
no longer a sharp observable. Rather, it is a thermal expectation value that mixes contributions from different states 
in the fermion Hilbert space. This is somewhat analogous to the recent observation of Jackiw et al |Q that fractional 
electron number in liquid helium bubbles is not sharp, due to state mixing. 

In the sigma model case, the resummed derivative expansion result (Eq) for (A^)t suggested a natural separation 
of the finite temperature induced fermion number (iV)T into a zero temperature topological piece (iV)o and a non- 
topological piece that represents the plasma response to the background field (in this case, an inhomogeneous electric 
field). This separation is, in fact, a general feature of the induced fermion number: (7V)t splits into a topological 
piece {N)q that arises due to vacuum polarization effects (and only involves the ground state of the second quantized 
Hilbert space), and a nontopological piece that can be interpreted in terms of the thermal occupation numbers of 



Paired states 




111 SHC B 



FIG. 8: The general structure of the Dirac spectrum for a kink background. There is a single unpaired bound state at _B = —m. 
(This assumes ^ > 0. If ^ < 0, then the unpaired bound state would be at i? = -\-m.) There are continuum thresholds at 

E = ±m secO — ±\/m'2 + 4>l. There may or may not be additional bound states with energies in the range m < |_E| < msec6, 
but if these states are present they are necessarily paired in ±E pairs. 
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-ni m 

J ^ 



FIG. 9; The structure of the Dirac spectrum for a sigma model background. The only general thing that can be said is that 
the continuum thresholds appear at i? = ±m. There is no particular symmetry to the spectrum, and there may or may not be 
one or more bound states in the gap. This depends sensitively on the detailed form of the angular field 0{x). 



nonvacuum states. The topological piece is temperature independent, while the nontopological piece is temperature 
dependent. To see how this separation arises, note the simple identity 

tanh (^^^ = 1 - 2n{E) (43) 

where the Fermi occupation number distribution is 

Then the finite temperature induced fermion number ( p^ separates as 

{N)t ^ -t: dEa{E)sign{E) + / dE a{E) sign{E) n{\E\) (45) 

^ J-oo J-oo 

The first term in (|4|) is just the zero temperature fermion number (||), and is topological. However, the fact that the 
spectral asymmetry integral in (^) is always topological does not necessarily mean that i7odd(-E) is itself topological. 
The second term in ( ^ ) is generically nontopological due to the weighting of the integral with the Fermi distribution 
factor. However, in the special case of a kink background, crodd(_E) is itself topological [see the Callias index theorem 
result (p^)]. Thus the temperature-dependent second term in (|4^) is also topological in the kink case, even though 
the integral involves the Fermi factor. 

This separation can be understood in terms of the Dirac spectrum of the fermions in the given background. For 
example, in the kink case, the Sommcrfeld- Watson integral representation (|l^) for the induced fermion number can 
be re-expressed using (|l9|), ( ^ and (^ ) as: 

, , 2m2tan^ f°° , n(E) , , 

(iV)T = (iV)o-sign(6l)n(m) + / dE ) ' = (46) 

{E"^ — m?)y E"^ — m^sec^O 



msccy 



For any kink background (with 6 positive), the spectrum of the Dirac Hamiltonian (ra) always has the form shown in 



Fig. 8. There are continuum states for \E\ > y + 4'2 — 'msecO, where we recall that 9 = arctan((/i2 /^ti) , and (j>2 is 
the asymptotic value of the kink field at a; = +oo. In addition, there is always a bound state ai E = — sign(0) to: 

tpti = ( r ^ fx I 1 ) with energy E = — m, if ^ > 



f 




i/,^ = ('^^^^^ "^^^ \ with energy E = +m, if ^ < (47) 



There may or may not be additional bound states with m < \E\ < m sec 0. If these additional bound states are present 
they necessarily occur in ±E pairs, because of the quantum mechanical SUSY of the Dirac Hamiltonian (||) in the 
kink case. Thus, the contributions of these paired bound states to the second integral in (|45| ) cancel in pairs, while 
the unpaired bound state at i? = — sign(^) to leads to the — sign(^)n(TO) term in (^). The remaining integral over 
the continuum states leads to the integral term in (^), as an integral of the Fermi factor n(|i?|) over the continuum 
beginning at the threshold energy toscc^ with the integrand involving the odd part of the spectral function derived 
from the Callias index theorem result ([l4|). 

In the sigma model case, the Dirac spectrum is very different. As discussed in Section III.B, the chiral field 9{x) is 
equivalent to a spatially inhomogeneous electric field E{x) = —^9"{x) for fermions of mass to. For a slowly varying 
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background, with 9' <C m, the leading order effect on the spectrum is that of a local chemical potential ii{x) = —^0'{x). 
Thus, we can approximate as 



(iV)T « (iV>o + / dx 



OO /'OO 



dk ( 1 



_^ 27r ye'3(v"='+™^+«'/2) + 1 g/3(Vfc2+m2-e'/2) _,_ y 



OC 



^ (iV)o-^^e-™/^y da;sinh(^^j , as T^O (48) 

which agrees precisely with the low temperature resummed derivative expansion result (|2^), and which agrees with 
the plasma linear response result as explained in Section III.B. 

From numerical analysis of the Dirac equation for a sigma model background, the Dirac spectrum has the form 
shown in Fig. 9. Note that, unlike the kink case spectrum depicted in Fig. 8, there is no particular symmetry of the 
spectrum. The only general thing we can say is that the continuum thresholds are at i? = ±m. There may or may 
not be one or more bound states for \E\ < m. The existence of (and the precise location of) bound states is highly 
sensitive to the actual shape of the 9{x). As there is no special symmetry in the spectrum, there is nothing general 
that can be said about the odd part of the spectral function. Nevertheless, the first term in (^5|), which is {N)o, 
yields a topological result since the integral can be written, using Levinson's theorem, in terms of the phase shift at 
infinity, which is topological |jl^, |l^ [l^. This application of Levinson's theorem doesn't work for the second term 
in ( p5| ) because of the Fermi factor. This is another way to understand why the finite temperature induced fermion 
number is generically nontopological. At low temperature, we see from (^ ) that the dominant correction to the zero 
temperature fermion number is determined by the bound state energy with the lowest magnitude: 

(Af)T-(iV)o+sign(^,ni„)e^''l^"""l+... , as T^O (49) 

Since the location of i^min is sensitive to the detailed shape of 9{x), this shows that (A^)t is generically nontopological. 
If there is no bound state, then the leading correction comes from the threshold and is ~ e~™^, with nontopolgical 
subleading corrections. 



The fluctuation (33) can also be expressed in terms of the Fermi distribution functions n{E) in ([14[). Noting the 
simple identity 

Isech'(^) =n{E){l~n{E)) (50) 



4 V 2 

we see that (IS^) can be written as 



{AN)^= dEa{E)n{E) {l-n{E)) (51) 

This formula is natural, since it is well known | |45[ | that for noninteracting fermions, the fluctuation in the occupation 
number of the state with energy E is {n\) — (nE)'^ = (it-e) (1 ~ ("■£;))■ In this paper we are considering the fermions 
to be in a fixed static background, so we are simply populating the single-particle energy levels of the corresponding 
Dirac Hamiltonian with noninteracting fermions, according to Fermi-Dirac statistics. This explains why the fluctuation 
vanishes at zero temperature, as the rt(l — n) factors vanish exponentially fast. It also explains why the fluctuation 
is nontopological, since the integral in (|l]) requires detailed knowledge of the spectrum. 

For the family of special kink backgrounds in (jsj), the integral representation result (^ ) for the fluctuation can be 
re-expressed in terms of the Fermi occupation numbers of the states in the spectrum: 

2m|tan^| E n{E) {I - n{E)) 



9rrj tan /9 

(A7V)2 = n{m){l~n{m)) + 2^n{Ei){l-n{Ei)) ^ ^ / dE 

1—1 ^ mscc^ 



4m|tan^| ^ / A E n{E) {\ - n{E)) ^^2) 

where the paired bound state energies zkEi for the special kink backgrounds (^) are at 

Ei=mJl + {l~^-^)tan''9 , / = 1, 2, . . . , (j - 1) (53) 
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This expression (|5^) can be interpreted directly in terms of the general expression (51) : we can easily identify the 
separate contributions from the unpaired bound state &t E — —to, from the paired bound states &t E = ±Ei, and 
the contributions from the continuum cuts. Since the integral in (^l]) involves the even part of the spectral function, 
the symmetry of the spectrum indicated in Fig. 8 does not lead to any particular simplification for the fluctuation, as 
there are no cancellations. Hence, (AiV)^ is nontopological, as it is manifestly dependent on the scale parameter j. 

In the sigma model case, we can interpret the equation (|5l|) as follows. As discussed previously, in the derivative 
expansion, 6' ^ to, limit, the leading effect of the background is that of a local chemical potential ii{x) — —^6'(x), in 
which case the local Fermi distribution factors are given by (p9|). Thus, the fluctuation can be computed from (O) as 



(AiV)l 



dx I — [/+(1 - /+) + /-(!-/-)- 2/o(l - /o) 

OO 

nl 



oo 



k^ + to^ — fJ-{x)) ^1 — n{\/ k^ + m? — ii{x))j 



+n{\/ k^ + + fJ,{x)) ( 1 — n{yk^-\-rn^ + fi{x))] — 2n{\/ k^ + m^) ( 1 — n{\/ k^ + mn?) 



(54) 



where the last term is the subtraction of the free case with /q — n{\/k'^ + m?). At low temperature, ( |54|) gives the 
leading contribution 



(AiV) 



dx 



dk 
2^ 



(55) 



where we have used the low temperature behavior of the k integral in (^l|). This expression ( ^5| ) agrees precisely with 
the resummed derivative expansion result (H3) for the fluctuation. 



VI. CONCLUSIONS 



To conclude, we reiterate that the finite temperature induced fermion number {N)^ is generically nontopological, 
and is not a sharp observable. The induced fermion number naturally splits into a topological temperature-independent 
piece that represents the effect of vacuum polarization on the Dirac sea, and a nontopological temperature-dependent 
piece that represents the thermal population of the available states in the fermion spectrum, weighted with the 
appropriate Fermi-Dirac factors. As the temperature approaches zero, the nontopological terms vanish exponentially 
fast, and we regain smoothly the familiar topological results at zero temperature. But at nonzero temperature, the 
nontopological contribution is more sensitive to the details of the spectrum, and so is generically nontopological. This 
is illustrated explicitly by a derivative expansion calculation, which is resummed to all orders at low temperature, 
for the case of a sigma model background. The induced fermion number (A^)t only depends on the odd part of 
the spectral function. This means the case of a kink background is special because the Dirac spectrum for a kink 
background is remarkably symmetric, and the odd part of the spectral function is itself a topological quantity, in the 
sense that it only depends on the asymptotic value of the kink field, rather than on the full details of the shape of 
the kink field. Thus, the kink case is non-generic, and the finite temperature induced fermion number is actually 
topological, albeit a complicated function of the temperature and the asymptotic value of the kink field. On the other 
hand, the fluctuation, (AA^)t, is determined by the even part of the spectral function, for which no general topological 
information is known. Thus, the fluctuations are always nontopological at finite temperature. Wc illustrated this 
explicitly with an exact solution for a special class of kink backgrounds, and with a derivative expansion calculation for 
the sigma model case. The fluctuation vanishes at zero temperature, which shows that the induced fermion number 
is a sharp eigenvalue at zero temperature, in agreement with the arguments of |2^, ^ |3^. But at finite 

temperature, the nonvanishing of the fluctuations indicates that the induced fermion number is not a sharp eigenvalue 
at nonzero temperature. The source of the fluctuation is clearly the dispersion introduced by mixing states in the 
thermal average. An analogous example of an induced fermion number expectation value that is not sharp has been 
discussed recently in Q. 

Some of these results generalize immediately to higher dimensional cases. For example, the dependence of (A^)t on 
just the odd part of the spectrum means that if the background has the same SUSY isospectral features as the kink 
case, then the finite temperature induced fermion number will be topological. This was in fact found to be true in 
where the 3-1-1 dimensional induced fermion number in the presence of a static SU{2) 't Hooft-Polyakov monopole was 
shown to have exactly the same form as in the 1-1-1 dimensional kink case, with the identification of the asymptotic 
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value of the kink with the asymptotic value of the magnitude of the monopole's Higgs field. This gives an interesting 
finite temperature remnant of the quantum mechanical SUSY present in the Dirac spectrum. Conversely, without 
such a special symmetry of the Dirac spectrum, the finite temperature induced fermion number will be generically 
nontopological. So, for example, for a chiral SU{2) sigma model with a Skyrme background in 3 + 1 dimensions, it is 
well known from numerical work that the Dirac energy spectrum of the fermions is not symmetric, and the energy of 
a possible bound state (and indeed the number of bound states) is highly sensitive to the details of the shape of the 
radial hedgehog field p9[ . Since the leading low temperature contribution to the finite temperature fermion number 
in (^) is determined by the lowest magnitude bound state energy [see (|49|)], this shows that in this case the finite 
temperature induced fermion number is nontopological. Likewise, the fact that the finite temperature fluctuation 
(AiV)x depends on the even part of the spectral function indicates that this will be generally nontopological. Finally, 
the interpretation in terms of the fluctuations for individual energy levels in the spectrum shows that the fluctuation 
is nonzero at finite temperature, but vanishing at zero temperature. Thus, the induced fermion number is a sharp 
observable at zero temperature, but not at finite temperature. In particular, for a Skyrme background in 3 + 1 
dimensions, the finite temperature fermion number is not simply the (topological) winding number of the Skyrme 
field, but a much more complicated nontopological object. 
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